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Variational View to Optimal Stopping 
Problems for Diffusion Processes and 
Threshold Strategies* 

V.I. Arkinj A.D. Slastnikov^ 


Abstract: We describe a variational approach to solving opti¬ 
mal stopping problems for diffusion processes, as an alternative 
to the traditional approach based on the solution of the free¬ 
boundary problem. We study smooth pasting conditions from 
a variational point of view, and give some examples when the 
solution to free-boundary problem is not the solution to optimal 
stopping problem. A special attention is paid to threshold strate¬ 
gies which allow reduce optimal stopping problem to more simple 
one-parametric optimization. Necessary and sufficient conditions 
for threshold structure of optimal stopping time are derived. We 
apply these results to both investment timing and optimal aban¬ 
don models. 


1 Introduction 

Let Xt, f > 0 be a diffusion process with values in D C defined on a 
stochastic basis (fl, A, {Xt, t > 0}, P). 

Let us consider an optimal stopping problem for this process: 

U{x) = sup (1) 
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where g : D —)■ is payoff function, p > 0 is discount rate, and E* 

means the expectation for the process Xt starting from the initial state x. 
The maximum in (1) takes over some class Xi of stopping times (s.t.)^ r. 
Usually (in classic case) Xi is the class of all s.t. with respect to the natural 
filtration Xf = a{Xs, 0 < s < t}, t >0). 

The above problem is one of the classical problems in stochastic control 
theory. It has a long and wide history in literature and applications. 

There are two main approaches to solving an optimal stopping problem 
for diffusion processes. 

The first one, usually referred as Markovian (or ‘mass’), embeds underly¬ 
ing optimal stopping problem into the family of problems (1) with all possible 
initial states x of the process Xt- In this case to solve problem (1) means to 
find the value function U{x) as well as stopping time t*{x), at which sup in 
(1) is attained (see, e.g., [20, 19]). 

It follows from general theory that U{x) is the smallest excessive (more 
precisely, p-excessive) majorant of payoff function g{x) (Dynkin’s characteri¬ 
zation, see [14]). Developing this approach in one-dimensional case, Dayanik 
and Karatzas [11] gave a characterization of the value function of the optimal 
stopping problem (1) as the smallest nonnegative majorant of payoff g{x), 
which is concave regarded the certain function which is specified by local 
characteristics of a given diffusion process. Presman [21] proposed how to 
derive value function using sequential modification of payoff function. 

Knowing the value function U{x) one can find the optimal stopping time 
T*{x) as the first exit time of process X^ out of the ‘continuation set’ C = 
{x : U{x) > g{x)}. 

The traditional method for hnding the value function for optimal stopping 
problem (1) is based on solving the related free-boundary problem; to find 
unknown function H{x), x E D and set U C R'^ such that 


L,H{x) = pH{x), X E C; 

(2) 

H{x) = g{x), X E dC; 

(3) 

grad 77(x) = gia.dg{x), x E dC 

(4) 


where L is the infinitesimal generator of Xt, dC is the boundary of the set C. 
The condition (3) is called “continuous pasting”, and (4) - “smooth pasting” 
condition. 

The solution to the above free-boundary problem is considered as a can¬ 
didate for solution to optimal stopping problem. And then, using some verifi¬ 
cation arguments it is necessary to prove that solution {H{x), C) to the free- 

^In this paper we consider stopping times which can take infinite values (with positive 
probability) 
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boundary problem really provides the solution of the optimal stopping prob¬ 
lem. Namely, H{x) will be a value function and t*{x) = inf{f > 0 : X* ^ C*} 
will be an optimal stopping time. This free-boundary approach is described 
in a lot of papers and textbooks, for example, in the monograph by Peskir 
and Shiryaev [20]. 

Another approach deals with solving an optimal stopping problem (1) 
for hxed initial state Xq = x. And, hrst of all, it should be noted the 
martingale approach by Beibel and Lerche [7, 8]. Their basic idea is to 
represent the functional in the right-hand of optimal stopping problem (1) 
as a product of positive martingale and ‘gain’ function. And then, using 
martingale arguments it can be shown that maximization of the gain function 
gives (under enough weak assumptions) a solution to initial optimal stopping 
problem. Other applications of martingale methods to solve optimal stopping 
problems one can hnd, e.g., in [20, Ch. Ij. 

The present paper develops so-called variational approach to solving an 
optimal stopping problem for hxed initial state of the process X*, described 
in [5]. In the framework of this approach we propose to hnd a solution to the 
problem (1) over the class Ai of stopping times, which are the hrst exit times 
of the process X^ from the sets belonging to the given family, and to make 
optimization over this family of sets. As an argument for such a reduction 
of the class of stopping times may be the fact that under enough general 
assumptions an optimal stopping time in problem (1) can be hnd as the hrst 
exit time of the process Xt out of the open set C = {f/ [x) > g{x)}. Hence we 
can take Ad as a family of hrst exit times from all open sets (in R'^). Since 
in one-dimensional case any open set can be represented as countable union 
of disjoint intervals, then optimal stopping problem can be reduced (in this 
case) to hnding optimal hrst exit time from intervals {a,b), I < a < b < r, 
which contains starting point x of the process Xt. Necessary conditions for 
optimality of such an interval were obtained in [Ij. The above mentioned 
intervals are an example of the fact that if the family of sets is chosen ‘well’, 
then maximization in hrst exit times from these sets can give a solution to 
optimal stopping problem over all stopping times. And conversely, one can 
formulate ‘an inverse problem’; Under hxed family of sets to describe all 
conditions on the process and the payoh function, under which an optimal 
solution (over all stopping times) is the hrst exit time of the process from a 
set belonging to the given family. 

Note, that for a lot of optimal stopping problem, especially in multi¬ 
dimensional case, it is impossible to derive an explicit solution, while the 
exact solution is very hard for calculations and not suitable for further anal¬ 
ysis. Thus, if hnding an optimal stopping decision is not the hnal goal of 
study (for example, in investment models, see [4]), then it makes sense to 
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restrict considerations to simple class of stopping times in order to obtain 
any ‘reasonable’ solution which will be tractable and suitable for analysis. 

At last, the variational approach gives a different look (compared with 
traditional one) to smooth-pasting principle (4). In the framework of this 
approach a smooth-pasting condition can be viewed as first-order optimality 
condition for a certain function, while an optimal stopping strategy is as¬ 
sociated with a maximization of that function. Thus, a difference between 
stationary point and point of maximum can give some non-exotic examples 
when there are many solutions to free-boundary problem and solution to free¬ 
boundary problem does not give a solution to optimal stopping problem. 

The paper is organized as follows. In Section 2 we describe a variational 
approach for solving an optimal stopping problem. For the case, when the 
underlying class of optimal stopping is first exit times from one-parametric 
family of sets in an optimal stopping problem can be reduced to one¬ 
dimensional maximization of some function. For this case we give necessary 
and sufficient conditions for optimality of stopping time over the class under 
consideration. 

Frequently, it is optimal to stop when the process exceeds some level 
(threshold strategy). Similar threshold decisions arise, e.g., in mathematical 
finance [22], investment models under uncertainty (real option theory) [13], 
etc. Almost all known decisions in real options theory have a threshold 
structure (see [13]). For example, solutions ‘to invest or not’, ‘to abandon or 
not’ depend on whether the observed values (which determine a decision) will 
be more or less than some level (threshold). In Section 3 it is demonstrated 
how a variational approach works for one-dimensional diffusion processes and 
two classes of one-parametric sets (/- and r-intervals), which are generated by 
threshold strategies. We give here necessary and sufficient conditions under 
which the optimal stopping time will have a threshold structure. At the 
end of this section we apply the obtained results to two fundamental models 
in real options theory: investment timing model and optimal abandonment 
model. 

In Section 4 we return to general variational approach for solving an 
optimal stopping problem and consider free-boundary problem (for threshold 
case) from a variational point of view. Namely, we give an example when 
a solution to free-boundary problem is not a solution to optimal stopping 
problem. Moreover, using second-order optimality conditions we prove some 
results about a relation between solutions to free-boundary problem and 
optimal stopping problem. 
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2 Variational approach to optimal stopping 
problem 

In this section we develop approach to solving an optimal stopping problem 
which we shall refer as a variational (see [5]). In the framework of this 
approach, one can dehne a priori a class M of stopping times which are the 
hrst exit time out of the set (from a given family of sets), and hnd optimal 
stopping time over this class. Besides, unlike the mass setting of an optimal 
stopping problem, we study the individual problem (1) for the given (hxed) 
initial state of the process Xq = x. 

Let Q = {G} be a given family of regions in M”, tg = tg(x) = inf{t > 
0 : Xt ^ G} be a hrst exit time of process Xf out of the region G (obviously, 
Tg = 0 whenever x ^ G), and Ai(G) = {tg, G G be a set of hrst exit 
times for all regions from the class Q. 

Under hxed initial value x for any region G G ^ we dehne the following 
function (of sets) 

Vg{x) = E^giXra)e-^^l^ra<oo}. (5) 

Outside the region G this function equals payoh function g, and inside G the 
function Vg{x) can be derived (under some weak assumptions) as a solution 
to boundary Dirichlet problem (see, for example, [14], [19]); 

L,u{x) = pu{x), X G G; , . 

u{y) g{x), y e G, y ^ X e dG. 

In order to calculate functions of the type (5) one can use martingale 
methods also (see, for example, [20], [12]). 

Thus, a solving an optimal stopping problem (1) over a class of s.t. Xi = 
X4{Q) can be converted to a solving the following variational problem; 

1 / 0 ( 0 ;)sup. (7) 

G&g 

If G* is an optimal region in (7), then the hrst exit time from this region 
Tg* will be the optimal stopping time for the problem (1) over the class 

M=M{g). 

If the class of regions Q is chosen ‘well’, it is possible to prove that s.t. 
Tg* will be also an optimal stopping time for problem (1) over all s.t. A4. In 
[5] such approach was realized for two-dimension geometric Brownian motion 
Xt and homogeneous payoh function g. 

A close approach is developed in [1], where an optimal stopping prob¬ 
lem for one-dimensional dihusion is solved by mathematical programming 
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technique. However those method uses a few properties of one-dimensional 
diffusion, which are not valid in multi-dimensional case. 

As for multi-dimensional processes, the hrst-order conditions as a heuris¬ 
tic method for hnding boundaries of optimal “continuation sector” in optimal 
stopping problem for bivariate geometric Brownian motion and homogeneous 
(of hrst order) payoff function was used in [12]. 

Let us note, that the calculation of the optimal stopping time over a 
given class of regions represents, to our opinion, a practical interest. Indeed, 
free-boundary problem for multi-dimensional diffusion processes has no (as a 
rule) explicit solution. Therefore, it has a sense to restrict our consideration 
to more simple stopping times and corresponding regions, for which it is 
possible to derive the function of sets Vg{x). Also, numerical methods can 
be applied for solving the problem (7) with hxed initial state Xq = x. 

An idea of variational approach is general enough and can be applied not 
only for a diffusion processes and payoff functions of the type (1). 

2.1 One-parametric family of regions 

Under some additional assumptions a general variational problem (7) can be 
simplihed and be made more convenient for study. 

Let Q = {Gp, p G P C be one-parametric family of regions in R”, 
Tp = inf{t>0 : Xt ^ Gp}, 

V{p-, x) = Vg,{x) = (8) 

(see formula (5)). 

The function V{p]x) is dehned on P x D, and, obviously, V{p;x) = g{x) 
for X ^ Gp. 

Further, we assume that a family of regions {Gp} satishes the following 
conditions: 

(Al) Monotonicity: Gp^ C Gpj whenever pi < P 2 - 

(A2) Thickness: Every point x E D belongs to the boundary of the unique 
set from the class Q. The parameter of those set will be refered as q{x), 
so X E dGq(x)- 

Under a thickness property one can write: 

V{q{x)-,x) = g{x) \/xeD. (9) 

For a one-parametric case the variational problem (7) can be rewritten 
as one-dimensional optimization: 

U(p; x) —)■ sup, (10) 

pGP 
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where V{p;x) is specified in (8). 

Under the assumptions (Al) and (A2) necessary and sufficient condi¬ 
tions for maximization of V (p; x) in p are given by the following 

Theorem 1. 

i) If p* =p*{x) is the solution to the problem (10) then conditions 

V (p; x) <V (p*; x) whenever p < p*, x e GpU dGp, (11) 

V (p; x) < g{x) whenever p > p*, x E Gp \ Gp* (12) 

hold. 

ii) If for some p* = p*{x) 

V{pi;x) > V{p 2 ,x) whenever p* < pi < P 2 , x G Gp^ (13) 

and condition (11) hold, then p* is the solution to the problem (10). 

Proof, i) Inequality in (11) is the direct consequence of optimality of 
p* in problem (10). For p > p* and x E Gp\ Gp* we have g{x) = V (p*; x) > 
U(p;a;), i.e. (12). _ 

ii) Show that V{p;x) < V{p*;x) for any p E P. Denote Gp = GpU dGp 
— the closure of the set Gp. 

Let p < p*. If a; ^ Gp*, then x ^ Gp (due to monotonicity of regions) 
and, therefore, V{p;x) = g{x) = V(p*;x). If a; G Gp* and x E Gp, then 
V{p-,x) < V{p*-,x) (due to (11)). Finally, if a; G Gp* and x ^ Gp, then 
q{x) < p*, X E dGq(^x)i and using (9), (11), we have: V{p-,x) = g{x) = 
V{q{x);x) < V{p*-,x). 

Consider the case p > p*. If a; ^ Gp, then x ^ Gp* (due to monotonicity of 
regions), hence, V{p;x) = g{x) = V(p*;x). Whenever x E Gp and x E Gp*, 
then V(p;x) < V(p*;x) due to (13). When x E Gp and x ^ Gp* one can see 
that p* < q{x) < p. Therefore, (13) implies: V(p;x) < V{q{x)-,x) = g{x) = 
V{p*;x). □ 

As one can see there is a gap between necessary conditions (11), (12) and 
sufficient conditions (11), (13). But for one-dimensional diffusion processes 
it will be derived that the their direct consequences give necessary and suffi¬ 
cient conditions for optimality of stopping time, so the mentioned above gap 
disappears (see Theorem 2 in Section 3). 

2.2 A variational look to smooth pasting principle 

The variational approach can give a new look to a smooth pasting principle, 
which is crucial in solving free-boundary problem. 
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Let the set of states D C M"- of the diffusion process Xt be an open set, the 
assumptions (A1)-(A2) on one-parametric family of sets Q = {Gp, p G P}, 
where P is an open set in hold. Suppose that functions g{x), as well as 
q{x) (the parameter of the region whose boundary passes through the point 
x), are differentiable in all arguments. Note, that the function q{x) will be 
smooth, for example, in the case when regions’ boundaries dGp are specihed 
by surfaces of the type {'l/(p, y) = 0, y E M”}, where is continuously 

differentiable in [p, y) and 'Lp(p, y) is non-zero. Moreover, assume there exists 
differentiable function F{p, x) on P x D such that F{p, x) = V (p; x) for 
p E P, X E Gp (where V{p;x) is dehned in (8)). 

Let p{x) be a stationary point of the function F{p,x) in p, i.e. 
Fp{p{x),x) = 0 (x G D). The continuous pasting condition V{q{x);x) = g{x) 
(see (9)) implies 

Fp{q{x),x) gradg(a;) + F^{q{x),x) = gradg{x), x E D. (14) 
Thus, if a; G dGp(^x)i then q{x) = p{x) and, therefore, 

F'xi.Pi.x),x) = giadg{x), x E dGp(^^). (15) 

This equality can be viewed as a variant of smooth pasting condition at 
the boundary of the set Gp(^x), whose parameter is a stationary point of a 
function F{p, x). 

Note, the set of such x that (15) holds can be empty. Consider the case, 
when stationary points p{x) = p do not depend on x. In this case the set 
of such X that relation (15) valid, is not empty. As we see below at Section 
3, such a situation emerges, in particular, for one-dimensional diffusions and 
the classes of /-intervals or r-intervals. In these cases the function F{p,x) 
has a multiplicative structure: F{p,x) = Fi{p)F 2 {x) — see (20), (26). 

Dehning the function F(x) = F{p,x), the relation (15) can be written as 
follows: 

gradh^(a;) = gradF(a;) = gradp(a;), x G dGp. (16) 

Taking into account that F{x) for a; G Gp is a solution to Dirichlet prob¬ 
lem (6), the equality (16) is a traditional smooth pasting condition, and, 
therefore, the pair {F{x), Gp) is a solution to free-boundary problem (2)-(4). 
Let us note, if gradg(a;) ^ 0 for the family of regions {Gp} and x E D, then 
(as one can see from (14)), the smooth pasting condition (16) is equivalent 
to stationarity of F{p, x) in p at the point p. 

On the other hand, if maximum of the function F{p, x) in p is attained 
at the point p* E P, then the hrst exit time Tp* from the region Gp* will be 
a candidate for an optimal stopping time over the class Ai(Q) (see Theorem 



Of course, point of maximum is a stationary point of the function F, but 
not vice versa. Hence, for such a case a solution to free-boundary problem 
can not give a solution to optimal stopping problem. We continue to discuss 
this question in Section 4. 

3 One-dimensional diffusion processes. 
Threshold strategies 

Let Xt be a diffusion process with values in the segment F C with bound¬ 
ary points I and r, where —oo < I < r < -|-cxd, open or closed (i.e. it may be 
(/,r), [l,r), or [l,r]), and its inhnitesimal generator has the following 

type : 

hf{x) = a{x)f'{x) + i a^{x)f''{x), (17) 

where a : F —)■ and a : F —)■ R),_ are the drift and diffusion functions, 

respectively. Denote X = (/,r). 

The process Xt is assumed to be regular; this means that, starting from 
an arbitrary point a; G X, this process reaches any point y eX in hnite time 
with positive probability. To guarantee the regularity of Xf we assume that 
the drift and diffusion functions are satisfy the following local integrability 
condition; ^ 

f —< oo for some e > 0, (18) 

Jx-e {y) 

at any x E X (see, e.g. [16]). 

It is known that under regularity conditions (18), on the interval X, there 
exist (unique up to constant positive multipliers) increasing and decreasing 
functions 'ip{x) and y:>{x) with absolutely continuous derivatives, which are 
the fundamental solutions to the ODE 

X,u{x) = pu{x) (19) 

almost sure (in Lebesque measure) on the interval X (see, e.g. [16, Chapter 5, 
Lemma 5.26]). Moreover, 0 < ipix), (p{x) < oc for x E X. Note, if functions 
a(x), a(x) are continuous, then ip, ip E C^(X). 

Below we introduce two natural one-parametric family of sets underlying 
a variational approach for one-dimensional diffusions. 

3.1 /-intervals 

As the first one-parametric family of sets we take intervals of the type Gp = 
{x E D ■. X < p}, p E which we call /-intervals (emphasizing that 
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left end is fixed and equal to left boundary point of the process values D). 
/-interval Gp is [l,p) or {l,p) in dependence on / G -D or not. Obviously, the 
class of /-intervals satishes conditions (A1)-(A2) in Section 2.1. 

Let us dehne = inf{/>0 : Xt ^ Gp} = inf{/>0 '■ Xt > p} — the hrst 
time when the process Xt leaves Gp. We will call as threshold stopping 
time (hrst exit time over threshold p). Let p G X} be a class of 

all such threshold stopping times. 

Along with the above stopping time let us dehne the hrst hitting time to 
threshold; Tp = inf{/>0 : Xt = p}, p G (/,r). 

Then for the above class of /-intervals the function V{p-,x), dehned in 
previous section, has the following representation; 

Lemma 1. If x,p eX, then 


V (p; x) 


g{p)'f’{x)/'ip{p), for X < p, 
g{x), for X > p, 


where 'ip{x) is an increasing solution to ODE (19). 
Proof. Due to known formulas (see, e.g., [15, 9]); 




'ip{x)/'ip{p), for X < p, 
(p{x)/(p{p), for X > p. 


( 20 ) 


Therefore, for x < p, obviously, = Tp and using the above formula we 
have; 

= Wg(XT,)e-'^' = g(p)E"’e-''^- = 

that proves lemma. □ 

For one-dimensional case we can enhance Theorem 1 and give necessary 
and sufficient conditions (which are the same) for optimality over the class 
of /-intervals (or threshold stopping times). Let us dehne 


h{p) = g{p)lfj{p)- (21) 

In order to exclude trivial cases (when optimal threshold stopping time 
does not exist) we assume that g{xQ) > 0 for some Xq G X. 

Theorem 2. Threshold stopping time is optimal in the problem (1) 
over the class for all x eT, if and only if the following conditions 

hold: 


h{p) < h{p*) whenever p < p*; (22) 

h{p) does not increase for p > p*. (23) 
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Proof. Using Lemma 1 it is easy to see that in our case the condition 
(11) is equivalent to (22), and (12) is equivalent to (23). Therefore, by 
Theorem 1, optimality of implies conditions (22)-(23). 

On the other hand, let (22)-(23) hold. Take arbitrary pi,P 2 such that 
p* < Pi < P 2 and X < p 2 . If a; < pi then V{pi^x) = h{pi)'ip{x) > h{p 2 )'ip{x) = 
V (p 2 ', x). And if Pi < x < p 2 then V (pp x) = g{x) = h{x)'ip{x) > h{p> 2 )'ip{x) = 
V{p 2 ', x). Thus, the condition (13) holds and we can apply Theorem 1 again. 
□ 

So, the optimal threshold p* is a point of maximum for the function h{p). 
This implies the necessity (under minor assumptions) of the ‘generalized’ 
smooth-pasting principle. 

Corollary 1. Let Tp*, p* G X, be the optimal stopping time in the problem 
(1) over the class and function g{x) has one-sided derivatives p'(p*±0) 
at the point p*. Then the function v{x) = supU(p;a;) has one-sided deriva- 

p&X 

tives at the point p* and the following inequalities hold: 

g'{p*+t)) = U(p*+0) < U(p*-0) < p'(p*-0). 


Proof. From (20) it follows 


v{x) = V{p*;x) 


h{p*)'ip{x), for X < p*, 
g{x), for X > p*. 


(24) 


Since p* is a point of maximum for the function h{p), then h'{p*—0) > 
0 > h'(p*-|-0). Therefore, 


v’{p*-0) = h{p*)i/j’{p*) = g'ip*-0) - h’{p*-0)i/j{p*) < g'{p* - 0), 

v'{p*+0) = g'{p*+0) = h\p*+0)ij{p*)+h{p*)ij’{p*) < h{p*)ij’{p*) = v’{p*-0). 
The corollary proved. 

Remark. If function g{x) is differentiable at the point p*, then the func¬ 
tion v{x) will be differentiable at the point p*, and v'{p*)=g'{p*). 

The necessity of the smooth pasting condition under some additional 
constraints on the process was shown in [20]. A result similar to ours was 
obtained in [23] 

The extended conditions of Theorem 2 will be necessary and sufficient 
for that a solution to stopping problem (1) over threshold stopping times 
remains optimal over all stopping times. 
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Let us define for p G X a restriction of a process Xt on = {x E D : 
X > p} as Yf = Xt/\Tp where Yq = x > p and Tp is the hrst hitting time of 
threshold p. In other words, the process Yt is started from a point greater 
than p and is absorbed when it reaches p. 

We will call the nonnegative function / ; —)■ M_|_ as p-excessive with 

respect to restriction of Xt on if for all a; > p and any stopping time r 

< fix). (25) 

where Tp = inf{t>0 : Xt = p}, p E {I, r). 

Theorem 3. Let Tp,, p* E X, he an optimal stopping time in the problem 
(1) over the class of threshold stopping times for all x eX, moreover, 
there exists g'{p*+0) and g{x) > 0 for x > p*. Then threshold stopping time 
Tp* is the optimal stopping time in problem (1) over class of all stopping times 
for all X eX, i.e. 


U{x) = snpWg{X^)e-P^l{r<oo} 

T 


h{p*)'ip{x), for X < p*, 
g{x), for X > p*, 


if and only if the function g{x) is p-excessive with respect to restriction of 
process Xt on D^*. 

Proof of Theorem 3. Let Tp* be optimal in problem (1) over all 
stopping times for all x E X. Then for x > p* and any stopping time r we 
have 

WgiX,)e-^%r<oo} < U{x) = g{x), 

where f = t A Tp*. It means that g{x) is p-excessive function with respect to 
restriction of process Xt on D^*. 

To establish the inverse implication let dehne the function 


<h(a;) = V{p*]x) 


h{p*)'ip{x), for X < p*, 
g{x), for X > p*. 


and show that E^<h(XT- < <h(a;) for any stopping time r and 

X eX. 

In order to prove it we use the following criteria of p-excessive functions 
from [17, Theorem 5.1], 

Lemma 2. A function F : D M+ satisfies the ineguality 


E^F{Xr)e-^l{r<oo} < F{x) 
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for all stopping times r and all initial states x ^ I, if and only if the following 
statements hold: 

i) F{x) is the difference of two convex functions for x G X; 

a) the measure CF{dx) = ^a‘^{x)F''{dx) + a{x)F'{x—0)dx — pF{x)dx, 
where F"{dx) means second distributional derivative, is non-positive; 

Hi) F is lower semi-continuous at absorbing boundary points. 

Apply this lemma to the function <h(a;). 

If a; < p* then <h(a;) = h{p*)if{x). Obviously, if{x) is the difference of 
two convex functions, because if'{x) has a hnite variation (at any segment 
[l',p*], I' > 1) and, therefore, is the difference of two increasing functions. 
Thus, condition i) of Lemma 2 is satished. Further, C^{dx) = h{p*)\hif{x) — 
pfj^xfdx = 0 a.s. (by dehnition of 'f{x)), i.e. ii) holds. At last, if / is 
absorbing state we can dehne if>{l) in continuity, so that condition iii) of 
Lemma 2 is also satished. 

If a; > p* then <h(a;) = g{x). Since g{x) is p-excessive with respect to 
restriction of process Xt on , then conditions i)-iii) of Lemma 2 are sat¬ 
ished. 

Finally, at the point x = p* we have £<h({p*}) = | o-^(p*)[p'(p*-T0) — 
h{p*)if'{p*)] < 0, since p* is the point of maximum for function h{p) (see 
Theorem 2) and, therefore h'{p*-\-0) < 0. 

Thus, all conditions i)-iii) of Lemma 2 are satished for the function 
<h(a;), X E D. 

Since is optimal over the class then Theorem 2 implies h{p*)if>{x) > 
h{x)if{x) = g{x) for x < p*, therefore <F(a;) > g{x) for all x. Using this in¬ 
equality and Lemma 2 we have 

E"p(X.)e-""l{.<oo} < E"$(X.)e-""l{.<oo} < <h(a;) 

for all stopping times r and x eX. Hence, U{x) < <F(a;). On the other hand, 
obviously, U{x) > <F(a;). 

Therefore, U{x) = <F(a;) = V(p*;x), i.e. Tp* is the optimal stopping time 
in problem (1) over all stopping times for all x. 

This completes the proof. 

As we see the question “Does optimal threshold stopping time Tp* remain 
optimal over all stopping times?” is reduced to problem of p-excessiveness of 
payoh function with respect to the restriction of underlying process on the 
range There are several conditions and criteria for p-excessiveness, see, 
e.g. [11, 17]. 

Now we give necessary and sufficient conditions for such a p-excessiveness 
of payoh function, which may be convenient for applications. 
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Statement. Let for a set of isolated points {ai,a ^,where 
p*<ai<a 2 <---<r, the funetion g'{p) he absolutely eontinuous on the intervals 
{p*, ai), (a*, Oj+i), i>l and there exist one-sided derivatives g{p*+t)), g\ai±tS), 
i>l, sueh that Y^a‘^{ai)\g\ai+tS)—g\ai—t))\ < oo. Then g{x) is p-exeessive 

i>l 

funetion with respeet to restrietion of proeess Xf on if and only if the 
following eonditions hold: 

i) hgi^p) < pg{p) a.s. (in Lebesque measure) for p > p*; 
a) g'{ai+0) — g'{ai—0) < 0 for all i >1. 

(This statement easy follows from Lemma 2.) 


3.2 r- intervals 

Now let us consider another one-parametric family of sets underlying the 
variational approach to solving optimal stopping problems, namely, the class 
of r-intervals of the type Gp = {x & D : x > p}, p &X. This class is similar 
to the class of /-intervals, but here we fix the right end of the range of the 
process D instead of its left end (as in /-intervals). 

In complete analogy with the above considerations we can dehne threshold 
stopping time rf = inf{/>0 : X* < p} as a first time when process Xt falls 
below threshold p, p G X, and the corresponding class 

Then for the class of r-intervals the function V{p-,x), dehned in (8), has 
the following representation (cf. (20)); 


V (p; x) 


g(p)(p(x)/(p(p), for X > p, 
g(x), for X < p, 


(26) 


where ip{x) is a decreasing solution to ODE (19). 

Thus, we can modify all results from the previous section for the class of 
r-intervals. 

Theorem 2'. Threshold stopping time is optimal in the problem (1) 
over the elass x eT, if and only if the following eonditions hold: 

g{p)/(p{p) < g{p*) / pip*) whenever p > p*; 
g{p)/p{p) does not decrease for p < p*. 


We can define = {x E D : x < p} for p E X and a restriction of a 
process Xt on as Yt = Xt/\Tp, where Yq = x < p and Tp is the first hitting 
time of threshold p. 
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In analogue with the case of /-intervals we call the function / : —)■ M+ 

/^-excessive with respect to restriction of on if for all a; < p and any 
stopping time r the inequality (25) hold. 

Now, the analogue of Theorem 3 is the following. 

Theorem 3'. Let p* G X, be an optimal stopping time in the problem 
(1) over the class threshold stopping times for all x E X, and there 

exists g'{p*—^)- Then threshold stopping time is the optimal stopping 
time in problem (1) over class of all stopping times for all x eX, i.e. 


U{x) = supE^ 5 ((X^)e ^^l{^<oo} 

r 


g{p*)^{x)/(p{p*), for X > p*, 
g{x), for X < p*, 


if and only if the function g{x) is p-excessive with respect to restriction of 
process Xt on X)p». 


Threshold stopping times rf appears, in particular, in optimal stopping 
problems with both integral and terminal payoffs: 




gi{Xt)e dt + go{Xr)e sup. 


(27) 


where go{x), gi{x) are given functions. 

Indeed, using strict Markov property of diffusion processes one can reduce 
the problem (27) to optimal stopping problem (1) with terminal payoff 


g{x) = go{x) — R{x), where i?(a;) = E"^ / gi{Xt)e dt. 


(28) 


By Green function (see, e.g., [9, 1]) one can get the following Green 
representation for R{x): 

R(x) = B-'(^ip{x) iiy)gi{y)H{y)dy + i’{x) J tp(y)g^(y)H(y)dyj , 

(29) 

where 5 = [f)'{x)(p{x)—ip{x)(p'{x)]/S'{x) > 0 is constant, H{y) = 2[a‘^{y)S'{y)] 
S'{x) = exp < — / ! dx > is derivative of the scale function for the pro- 

l J } 

cess Xt. 

If functions gi{x) and a{x), o'{x) are continuous, then Green representa¬ 
tion (29) implies that R{x) will be twice differentiable at interval X and 

R'{x) = B~^[ip'{x)Ii{x) + 'ip'{x)l 2 {x)], (30) 

R"{x) = B~'^[(p'{x)h{x) + fj'{x)l 2 {x)] - 2gi{x)/a‘^{x)-, (31) 


h{x)= ^lJ{y)gl{y)H{y)dy, hix) = (p{y)gi{y)H{y) dy, (32) 
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where B is defined in (29). 

Hence, Theorems 2', 3' together with the variant of Statement (for r- 
intervals) imply the following result concerning a threshold structure of so¬ 
lution to optimal stopping problem with both integral and terminal payoffs 
(27). 

Theorem 4. Let for some p* eX functions a{x), a{x), gi{x) be continu¬ 
ous, go{x) be twice differentiable and go{x) > R{x) on segment {l,p*]. Then 
threshold stopping time rf* = inf{t>0 : Xt < p*} is the optimal stopping 
time in problem (27) over the class of all stopping times (for all x E T) if 
and only if the following conditions hold: 

[go{p) - R{p)](p{p*) < [go{p*) - R{p*)](p{p) for p > p*- (33) 

h{p*)S'{p*) = g'o{p*)p{p*) - go{p*W{p*)] (34) 

Lfl'o(p) - P9o{p) < ILi?(p) - pR{p) for p < p*, (35) 

where function is specified in (32) and p{x) is a decreasing solution to 
ODE (19). 

Results presented in Sections 3.1 and 3.2 extend the corresponding results 
obtained in [2, 6, 10, 23]. 

3.3 Application to Real Options 

Investment timing problem. One of the fundamental problems in 
real options theory is to derive the optimal time at which an investor (or de¬ 
cision maker) should finance and launch an investment project — investment 
timing problem. 

Let / be a cost of investment required for beginning a project, and Xt is 
considered as Present Value (PV) from the project started at time t. As usual 
investment supposed to be instantaneous and irreversible, and the project — 
infinitely-lived. 

At any time a decision-maker (investor) can either accept the project 
and proceed with the investment or delay the decision until he obtains new 
information regarding its environment (prices of the product and resources, 
demand etc.). Thus, the main goal of a decision-maker in this situation is to 
find, using the available information, an optimal time for investing the project 
(investment timing problem), which maximizes the Net Present Value from 
the project: 

E {Xr - I) e-P" -E max, (36) 

r 

where the maximum is considered over all stopping times teM. 
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The majority of results on this problem which was introduced in [18], has 
a threshold structure: to invest when PV from the project exceeds the certain 
level (threshold). In the heuristic level it was shown for the cases when PV 
is modelled by geometric Brownian motion, arithmetic Brownian motion, 
mean-reverting process and some other (see [13]). And the general question 
arises: For what processes of PV from the project an optimal decision to an 
investment timing problem will have a threshold structure? Some sufficient 
conditions in this direction was obtained in [1]. We can give the necessary 
and sufficient conditions for optimality of threshold strategy in investment 
timing problem (36) as the consequence of our above results (but not a direct 
corollary) for the linear payoff function g{x) = x — I. We consider the case 
I < r else the optimal time in (36) will be -|-cxd. 

Theorem 5. Threshold stopping time Tp*, p* G (/,r), is optimal in the 
investment timing problem (36) if and only if the following eonditions hold: 

{p — I)fj{p*) < {p* — I)fj{p) iorp<p*; (37) 

^(p*) = (p*-J)^'(p*); (38) 

a{p) < p{p — I) forp>p*, (39) 


where if{x) is an inereasing solution to ODE (19) and a{p) is the drift fune- 
tion of the proeess Xt. 

Proof. Let r^*, p* G (/, r), is optimal in the investment timing problem. 
Then it is optimal in (36) over threshold stopping times This implies 

(due to Theorem 2) (37)-(38) and, besides, p* > I. According to Theorem 3 
g{x) = X — I is p-excessive function with respect to restriction of Xt on the 
range Up*. This imply (39) — see condition i) from the Statement. 

Now, let (37)-(39) hold. As in the proof of Theorem 3 consider the 
function 


<F(a;) = V{p*]x) 


{p* — I)'ip{x)/'ip{p*), foTx<p*, 
X — I, for a; > p*. 


Note, (37) imply <F(a;) > a; — / for all a; G X. 

It is easy to check (similar to the corresponding arguments in the proof 
of Theorem 3) that all conditions i)-iii) of Lemma 2 are satished. There¬ 
fore, E^$(X,)e-^"l|,<oo} < <h(a;) for all r and a; G X. Thus, E^(X, - 
I)e~^l{r<oo} < ^(a;) = V{p*]x), i.e. r^* is optimal in the problem (36) over 
all stopping times. 

Remark. As one ean see from the proof of Theorem 5, eonditions (37)- 
(39) are necessary and sufficient for optimality of stopping time Tp* in optimal 
stopping problem with ‘optional’ payoff function g~^{x) = {x — I)~^. 
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Let Xt be geometric Brownian motion with rate of drift a and volatility 
a. In this case, where {3 is the positive root of the equation 

— 1) + tt/S = p. If p > a then (3 > 1. The optimal threshold 

is p* = —— -I (see, e.g. [13, Ch. 5]). Note, if a ^ cxd then (3^1 and, 

therefore, p* —)■ cxd. Hence, a large volatility increases threshold for investing. 


Optimal exit (abandonment) problem Another fundamental prob¬ 
lem in Real options theory concerns hnding the optimal time at which a 
decision-maker, who receives a payoff from an acting hrm, should terminate 
operating and abandon hrm. 

Let Xt be the price (at time t) of the good produced by the hrm, and 
function gi{x) describes a dependence of hrm’s revenue on current price x. 

The optimal abandonment problem is, using the available information 
about current production prices Xt, to hnd a moment r for terminating of 
hrm’s activity (exit moment) such that the net present value of the hrm be 
maximal; 


yj dt — Le —)■ max, (40) 

where L > 0 is the abandonment cost. 

The above Theorem 4 allows to rewrite the necessary and sufficient con¬ 
ditions for optimality of threshold stopping time = inf{t>0 '■ Xt < p*} 
(with assumption L + R{x) < 0 for x < p*) in the abandonment problem as 
follows: 


[L + R{p)]lp{p*) >[L + R{p*)]lp{p) for p > p*; (41) 

h{p*)S'{p*) = W{p*y, (42) 

LR(p) > pR{p) + pL for p < p*, (43) 

where functions R{x) and hix) are specihed in (29) and (32) respectevely. 

Let us note that S'{p*) > 0, p>'{p*) < 0. Therefore, if follows from (42) 
that 

h{p*) = [ <p{y)gi{y)H{y)dy <0. (44) 

Jp* 

If we assume that function gi{x) is increasing, i.e. hrm’s revenue increases 
when current price rises, then from (44) we have immediately gi{p*) < 0. It 
means that a hrm should terminate an activity only when its revenue will be 
strictly negative. In particular, if gi{x) = x — c, where x is current price of 
production, and c is production cost, then optimal threshold of price p* for 
abandon a hrm should be less then production cost: p* < c. 
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If a firm under consideration operates, e.g., in oil exploration, the op¬ 
timal threshold p* can be interpreted as optimal (in sense of problem (40)) 
breakeven oil price, that characterizes how low price can fall before oil projects 
start shutting down. 

Consider the example (see [13, Ch. 7]), where diffusion Xf is geometric 
Brownian motion with rate of drift a < 0 and volatility a, and gi {x) = x — c. 
In this case, R{x) = x/{p — a) — c/p, 'p{x) = where /5i is the negative 
root of the equation — 1) + a/3 = p. 


/Si 


P 


a 


(c — pL). It is 


If c > pL, the optimal threshold is p* = 

/3i-l p 

easy to see that /3i > p/a, p* < c—pL and all conditions (41)-(43) hold. 
Note, that if volatility a tends to -|-oo, then /5i —)■ 0 and, therefore, p* —)■ 0. 
It means that large volatility implies low level of price for shut down and 
more long period before abandonment of hrm’s activity, even though current 
prohts Pt — c will be negative. 


4 Pree-boundary problem and 
threshold strategies 

Let us return to a variational approach, described at Section 2, and give a 
new look to a free-boundary problem for the ‘threshold case’. Consider one¬ 
dimensional diffusion process Xt with values in interval D = {l,r) and the 
class of /-intervals. Assume that both drift and diffusion functions of X^ are 
continuous, and payoff function is differentiable (on D). 

For this case a free-boundary problem can be written as follows: to hnd 
threshold p, I < p < r and bounded twice differentiable function H{x), I < 
X < p, such that 


L,H{x) = pH{x), 1 < X < p; 

(45) 

H{p-0) = g{p). 

(46) 

1 

o 

(47) 

Conditions (45)-(46) hold for the function 


H{x) = h/p/ip^x), 1 < X < p 

(48) 

where i/i{x) is an increasing solution to ODE (19) and h{p) 

= 9(p)/V’(p). 


Thus, smooth pasting condition (47) at the point p is equivalent to g{p)'ip'{p) = 
g'{p)'ip{p), or h'{p) = 0, i.e. p is a stationary point of the function h. 
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On the other hand, as the results of Section 3 show, the optimal threshold 
must be a point of maximum of the function h. Of course, point of maximum 
is a stationary point of the same function, but not vice versa. 

In economic literature (e.g. in Real Options theory) it is common opin¬ 
ion that solution to free-boundary problem is always a solution to optimal 
stopping problem. 

To demonstrate a difference between solutions to optimal stopping prob¬ 
lem and free-boundary problem, consider the following example. 

Example. A solution to free-houndary problem ean not give a solution 
to optimal stopping problem. 

Let us consider geometric Brownian motion Xf = ^expltcj} (where Wt 
be standard Wiener process), payoff function g{x)={x — 1)^ -f x^ for a; > 0, 
and discount rate p = h^/2 [6 > 0). The function g is smooth and increasing 
for all h > 0. 

For this case the free-boundary problem is the following one; 

{ ^x‘^H"{x) \xH\x) = pH{x), 0 < X < p* 

H{p*) = g{p*) (49) 

H'ip*) = g'{p*) 

Remind the designations Vp{x) = 'EAg{Xrj,)e~^'^^, and Tp = min{t > 0 : 

Xt>p]. 

For (5 < 3 the free-boundary problem (49) has the unique solution; 
H[x) = Vi{x) = x^, p* = 1. However, the stopping time Ti (a hrst exit 
time over the level 1) is not optimal since Vp{x) = ((p — l)^/p'^ -|- l)a;'^ —)■ oo 
when p ^ oo (when S < 3), and Vp(x) f V{x) = 2a;^ (when 5 = 3) for any 
a; > 0. 

For 5 > 3 the free-boundary problem (49) has two solutions; 

H{x)=Vi{x)=x^, p* = 1 and (50) 

H{x)=Vpg{x)=h{ps)x\ p*=ps=6/{6-3). (51) 

Note that Vpg{x) > Fi(a;). Thus, one of the solutions to free-boundary 
problem ((50), corresponding to the boundary p* = 1) does not give a solu¬ 
tion to the optimal stopping problem (which there exists, in contrast to the 
previous case). 

The above remarks together with classical results on extremal problems 
(especially, second-order optimality conditions) allow to derive new results 
about the relation between solutions to free-boundary problem (45)-(47) and 
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optimal stopping problem (1) over the class M.=Ai\Yi of threshold stopping 
times generated by all /-intervals. 

If function g{x) be twice differentiable at the point p, then representa¬ 
tion (48) implies that = h{p)'ip"{p). On the other hand, g"{p) = 

h''{p)'ip{p) + h{p)%lj''{p) since h'{p) = 0. Hence, we have 

h'\p)%l){p) = g'\p) - H'\p-Q). (52) 

Similarly, if h^^\p) = 0, k = 1,... ,n — 1 for some n>2 then 

h^^\p)Pj{p) = g^^\p) - H^^\p-Q). (53) 

Due to Theorem 2 if Tp* is optimal in the problem (1) over the class 
then function h{p) has maximum at the point p*. Applying second-order 
condition for extremum and relation (52) we obtain the following 

Proposition 1. If Tp* is optimal in the problem (1) over the class 
M=M[^ and function g{x) is twice differentiable at the point p*, then the 
pair {H{x),p*), where H{x) = h{p*)f{x), is the solution to free-boundary 
problem (f5)-(f7) and H"{p*—t]) > g"{p*). 

The inverse relation between solutions can be state as follows. 

Proposition 2. Let a pair {H{x),p*) be a solution to free-boundary prob¬ 
lem (45)-(47), such that H"{p*—0)>g"{p*) and one of the following condi¬ 
tions holds: 

1) {H{x),p*) is the unique solution to free-boundary problem (45)-(47); 

2) if there exists another solution {H{x),p) to free-boundary problem 
(45)-(47), then H{x) > g{x) for I < x < p*, and either 

(2a) p < p*; or 

(2b) p > p*, H^^\p—t))=g^^\p), k = 2,...,n — 1, H^'^\p—0)>g^"'\p) 
for some n>2. 

Then Tp* is optimal in the problem (1) over the class 

Proof. From the dehnition ofp*, inequality H”[p*—t))>g”[p*) and (52) 
it follows that h'{p*) = 0 and h'{p) strictly decreases at some neighborhood 
of p*. 

Let condition 1) holds. It is easy to see that h'{p) > 0 for p < p* and 
h'{p) < 0 for p > p* (else h'{q) = 0 for some q p* and {h{q)f{x),q) is 
another solution to free-boundary problem (45)-(47)). Then application of 
Theorem 2 gives the optimality of stopping time Tp*. 

Now, let 2) holds. Then representation (48) immediately imply that 
h{x) = g{x)/f{x) < H{x)/f{x) = h{p*) for I < x < p*. 
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Let us prove that h'{p) < 0 for all p > p*. Indeed, if h'{pi) > 0 for some 
Pi > p*, then there exists p* < po < pi such that h'{po) = 0 and h'{p) > 0 
for all Po < p < pi- Therefore, {h{po)'ip{x),po) is another solution to free¬ 
boundary problem (45)-(47), and due to 2) h^^\po) = 0, /c = 2,...,n — 1, 
^*'”Hpo)< 0 for some n>2, that contradicts to positivity of h'{p) for po < p < 
Pi- 

Hence, h'{p) < 0 for all p > p* and due to Theorem 2 Tp* is optimal in 
the problem (1) over the class of /-threshold stopping times. □ 

Let us demonstrate, how Proposition 2 works for the example considered 
above. 

In the case <5 < 3, we can’t apply Proposition 2 since H"{1)=5{5 — 

In the case S > 3, for the solution (51) with the boundary point p* = ps 
(see (51)) we have 

H"(p,) = S(S - > S(S - l)pj-^ + ^ = 9"M- 

Besides, h’{x) — {x— — 3) (^ ^ '’ ^ ~ < 0 for <5 > 3 and x < p,, 

therefore, H{x) = h{ps)x^ > h{x)x^ = g{x) for x < ps. Since for the 
another solution (50) boundary point p* = 1 < ps, then the condition 2) of 
Proposition 2 holds. 

Hence, Tp^ is optimal in the problem (1) over the class of threshold stop¬ 
ping times 

Moreover, it can be shown (using Theorem 3), that Tpg will be optimal 
stopping time over the class of all stopping times also. Indeed, by direct 
calculations one can obtain that 

lLg{x) - pg{x) < -1.5-—- <0 for a; > ps, 

0 — 3 

and we can apply Theorem 3. 


References 

[1] Alvarez L.H.R. Reward functionals, salvage values, and optimal stop¬ 
ping. — Mathematical Methods of Operations Research, 2001, v. 54, p. 
315-337. 

[2] Arkin V.I. Threshold Strategies in Optimal Stopping Problem for One- 
Dimensional Diffusion Processes. — Theory of Probability and its Ap¬ 
plications, 2015, V. 59, No. 2, p. 311-319. 


22 



[3] Arkin V.I., Slastnikov A.D. Optimal time to invest under tax exemp¬ 
tions. — In; The Shyraev Festschrift. From Stochastic Calculus to Math¬ 
ematical Finance. Springer, 2006, p. 17-32. 

[4] Arkin V.I., Slastnikov A.D. The Effect of Depreciation Allowances on 
the Timing of Investment and Government Tax Revenue. — Annals of 
Operations Research, 2007, v. 151, No. 1, 307-323. 

[5] Arkin V.I., Slastnikov A.D. A variational approach to an optimal stop¬ 
ping problems for diffusion processes. — Theory of Probability and its 
Applications, 2009, v. 53, No. 3, p. 467-480. 

[6] Arkin V.I., Slastnikov A.D. Threshold Stopping Rules for Diffusion Pro¬ 
cesses and Stefan’s Problem. — Doklady Mathematics, 86 (2), 626-629, 
2012 . 

[7] Beibel M., Lerche H. R. A new look at optimal stopping problems 
related to mathematical hnance. — Statistica Sinica, 1997, 7, 93-108. 

[8] Beibel M., Lerche H. R. A note on optimal stopping of regular diffusions 
under random discounting. — Theory of Probability and Its Applica¬ 
tions, 2001, 45 (4), 657-669. 

[9] Borodin A., Salminen P. Handbook of Brownian motion - Facts and 
formulae. Birkhauser Verlag, Basel-Boston-Berlin, 2002. 

[10] Crocce F., Mordecki E.. Explicit solutions in one-sided optimal stopping 
problems for one-dimensional diffusions. — Stochastics, 86 (3), 491-509, 
2014. 

[11] Dayanik S., Karatzas L. On the Optimal Stopping Problem for 
One-Dimensional Diffusions. — Stochastic Processes and Their Appli¬ 
cations, 2003, V. 107, p. 173-212. 

[12] Gerber H.U., Shin E.S.W. Martingale Approach to Pricing Perpetual 
American Options on Two Stocks. — Mathematical Finance, 1996, v. 
6, p. 303-322. 

[13] Dixit A.K., Pindyck R.S. Investment under Uncertainty. Princeton; 
Princeton University Press, 1994. 

[14] Dynkin E.B. Markov processes. I. Springer, Berlin, 1965. 

[15] Ito K., McKean H. Diffusion Processes and Their Sample Paths. 
Springer-Verlag, Berlin, 1974. 


23 



[16] Karatzas L, Shreve S. E. Brownian motion and stochastic calculus. 
Springer-Verlag, 1991. 

[17] Lamberton D., Zervos M. On the optimal stopping of a one-dimensional 
diffusion. — Electronic Journal of Probability, 2013, 18, no. 34, p. 1-49. 

[18] McDonald R., Siegel D. The value of waiting to invest. — Quarterly 
Journal of Economics, 1986, v. 101, p. 707-727. 

[19] Oksendal B. Stochastic differential equations. Springer-Verlag, Berlin, 
1998. 

[20] Peskir G., Shiryaev A. Optimal stopping and free-boundary problems. 
Birkhauser, 2006. 

[21] Presman E. Solution of the Optimal Stopping Problem for One- 
Dimensional Diffusion Based on a Modification of the Payoff Function. 
— In; A.N. Shiryaev et al. (eds.), Prokhorov and Contemporary Prob¬ 
ability Theory, Springer Proceedings in Mathematics & Statistics 33, 
2013, p. 371-403. Berlin, Heidelberg: Springer-Verlag. 

[22] Shiryaev A. N. Essentials of Stochastic Finance. Facts, Models, Theory. 
World Scientific, 1999. 

[23] Villeneuve S. On the threshold strategies and smooth-fit principle for 
optimal stopping problems. — Journal of Applied Probability, 44 (1), 
181-198, 2007. 


24 



